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An analytical High Value Target (HVT) acquisition model is 
developed for a generic anti-ship cruise missile system. 
The target set is represented as a single HVT within a field 
of escorts. The HVT's location is described by a bivariate 
normal probability distribution. The escorts are repre- 
sented by a spatially homogeneous Poisson random field 
surrounding the HVT. Model output consists of the prob- 
ability that at least one missile of a salvo acquires the 
HVT, conditioned on the number of missiles in the salvo 
which penetrate the HVT area defense. In addition, the fall 
of multiple penetrators is modeled using a conditional 
multinomial probability distribution. The model's equations 
are used to solve for an optimal missile seeker range gate, 
given a probability distribution describing the location of 
the HVT within the targeted formation at the time the 
missile commences its search. Included in an appendix is a 
time-dependent model describing HVT location which provides 
for HVT movement during missile time of flight up to the 
moment of active search. 
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I. INTRODUCTION 



With, the advent of modern, sophisticated anti-ship 
cruise missiles (ASCM), opposing naval battle forces have an 
over-the-horizon (OTH) strike potential which, when actual- 
ized by appropriate tactics, will herald a new era in naval 
combat. As an opening remark the above may sound dated, for 
the cruise missile has arrived and the technology is aboard 
many combat fleets of the world. Nevertheless, the implica- 
tion of next-generation microelectronics and artificial 
intelligence aboard the ASCM of tomorrow presents the 
forward-looking military operations research analyst with 
promising opportunities for tactical improvement. 
Revolutionary advances in ASCM technologies will provide 
future systems with multi-target discrimination ability, 
reduced soft-kill susceptibility, and other significant 
improvements. Hence, tactical options should be evaluated 
well in advance of a new missile's arrival in the Fleet. 

In an attempt to destroy an opposing, -escorted high 
value target (HVT), the decision-maker possessing an ASCM 
OTH capability must choose when to strike, using appropriate 
force ( e. g. , total salvo size), a decision surely influenced 
by targeting accuracy, intelligence regarding opposing force 
composition, and first-strike criticality. Quantitative 
answers to the above decision problem are by nature probabi- 
listic ones. The decision-maker can be assisted by an asso- 
ciated estimate of HVT kill probability from a model. 

The multi-target scenario has been modeled quite exten- 
sively using elaborate Monte-Carlo computer simulation tech- 
niques which, even after considerable programming input 
effort and program execution time, yield estimates that 
often neglect the essential vagaries of the war-at-sea envi- 
ronment. Simulation has contributed much insight to the 
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decision problem at hand/ but is generally inappropriate as 
a tactical decision-aid because of its inflexibility. The 
so-called Bernoulli Trials model is a much simpler, purely 
analytical probability model which has been employed as a 
tactical decision-aid, [Ref. 1]. 

In the Bernoulli Trials model, HVT defense is essen- 
tially limited by the number of available fire-control chan- 
nels, each of which is awarded a deterministic or 
probabilistic number of ASCM kills. Targeting information 
is assumed to be accurate and all ASCM that penetrate the 
HVT area defense have equal probability of acquiring any of 
the multiple targets present. In an extension to the model, 
HVT acquisition by the ASCMs is biased upward, thereby 
endowing the missiles with multi-target discrimination 
ability. The model's equations are solved iteratively for 
an estimate of salvo size needed to inflict a firepower-kill 
on the escorted HVT with a specified probability of success. 
The driving parameters are the number of escorts, the effec- 
tiveness of HVT defense, and the required number of hits on 
the HVT to achieve a firepower-kill. 

The Bernoulli Trials model and its simulation counter- 
parts make simplistic assumptions regarding ASCM target 
acquisition. An accurate target acquisition model is an 
essential part of any model that attempts to provide the 
strike planner with a cost-effective estimate of requisite 
salvo size, where "cost" refers to the number of ASCMs 
employed. (While it is true that the prudent tactician may 
wish to launch "extra" missiles as a hedge, it is also 
likely that missiles will be in short supply; hence, he 
needs to know what basic number he is supplementing with 
"extras". ) This paper offers a new, purely analytical HVT 
acquisition probability model designed to enhance existing 
multi-target ASCM OTH models. It retains much of the 
simplicity of the Bernoulli Trials model, but: 
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• Reflects statistical fluctuations in targeting accuracy 
in a better way. 



• Incorporates the attacker's missile search area 
guides his range gating decisions. 



and 
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I I . THE MODEL 



A. DESCRIPTION 

1. Overview and General Assumptions 

The ANALYTICAL HVT ACQUISITION MODEL estimates the 
probability that at least one penetrator (an ASCM that 
survives the HVT area defense) acquires a single HVT within 
an escort field of low value targets (LVT). Additionally, 
the expected number of penetrators that acquire the HVT and 
her escorts ( ie. , the fall of shot) is estimated. The 
model's outputs serve as reasonable Measures of 
Effectiveness of the decision-maker's choices under uncer- 
tainty; namely, given the current tactical picture and 
opportunity to strike at the enemy, how many missiles should 
he employ and how are they best programmed to achieve a 
desired probability of HVT acquisition? 

The choices of launch bearing and ASCM search 
programming are influenced by the accuracy and time late of 
the targeting information and the missile time of flight to 
the target area. Where and when opportunities exist, anal- 
ysis of the offense's information state concerning the 
targeted formation may be profitable even when information 
is scarce and uncertain. The model assumes the following 
information is available to the strike planner: 

• A bivariate normal probability density function lor 
other appropriate probability density) describing the 
location of the HVT within an array of escorts. 

• An estimate of LVT density about the HVT. 



• An appropriate set of missile seeker sensor Sweep 
Widths. 



• An estimate of target acquisition probability, given 
target detection. 
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• A simplified estimate of ASCM reliability during the 
search, detection, and acquisition phases of the attack. 

• A velocity- time random vector of HVT/escort field motion 
( optional ) . 

The missile firing strategy used in the model is 
that of near- simultaneous firings (a short firing interval 
with respect to appreciable target movement) on a common 
launch bearing calculated to intercept the mean of the HVT 
location density. If the model user has a postulated motion 
vector of the enemy formation, the HVT density can be 
updated to incorporate the time delay from the generation of 
the previous targeting ellipse to the time of arrival of the 
salvo at the search area, in which case the salvo launch 
bearing should be calculated to intercept the mean of the 
updated HVT density at the time of missile search. 
Adjusting the HVT density and launch bearing becomes neces- 
sary when missile time of flight and enemy motion result in 
significant displacement of the enemy formation from its 
location at the time of launch. 

2. Modeling the ASCM System 

A generic ASCM system is modeled. Missile flight is 
straight down the launch bearing through the programmed 
search area. The assumption of zero missile navigation 
error is made since incorporation of the relatively small 
errors associated with a highly accurate system is not apt 
to heavily influence model results. (Deference is made 

however, to the exacting navigational accuracy required in 
especially long range firings. ) The missile seeker sensing 
device is essentially off ( ie. , receives no information), 
prior to ASCM arrival at the search area. An important 
assumption concerning the offensive firing is the notion of 
conditional independence between shots; each missile within 
a salvo is given the same launch bearing and search program, 
but retains functional independence ( autonomy) with regard 
to target acquisition. 
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The effective ocean area searched out by the ASCM is 
a function of operator input (range gate) and the missile 
seeker sensor Sweep Width(s) given by a Definite Range Law 
of Detection. The resulting area is a rectangle (or set of 
rectangles when multiple Sweep Widths are defined). In the 
existing detection environment, individual targets may be 
grouped by various homogeneous characteristics ( e. g. , size, 
active or passive signatures, shape, motion, etc. ). For the 
purpose of detection, similar targets are subject to a 
common sensor Sweep Width. The model may be used in the 
situation where the HVT possesses a detectable feature that 
distinguishes it from the LVT class. For example, a 
sensor's Sweep Width for an aircraft-carrier may be wider 
than its Sweep Width for a destroyer-escort because of the 
carrier's relatively large size. 

A target is detected if it lies within its associ- 
ated rectangle of effective search and is "swept" by the 
missile seeker sensor ( ie. , is in the sensor's field of view 
and within half a Sweep Width of the search axis). The 
model assumes the flightpath and the search axis of each 
missile are the same. The model's generic ASCM accomplishes 
its search in a matter of seconds, making search analogous 
to taking a snapshot photograph ( an aggregated detection 
process). Consequently, target motion is treated as negli- 
gible throughout successive searches by individual ASCM; a 
reasonable assumption if the interarrival times between 
missiles within the salvo are small. The assumption is made 
that each of the reliable ASCM in the salvo detect the same 
set of targets. This is because of their identical search 
programming, the Definite Range Law of Detection, and a 
short firing interval. 

Once a target or number of targets have been 
detected, the missile enters an acquisition phase in which a 
decision to acquire a particular target is made. The model 
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associates a conditional probability of acquisition, an 
input, with each target ( HVT and one or more LVTs) detected 
by the ASCM. Implicit in the detection and acquisition 
modeling is the assumption that, for purposes of acquisi- 
tion, detected targets are processed independently and in 
sequence from near to far within the missile seeker sensor's 
range gate. Therefore, acquisition of a particular target 
is influenced by two factors following detection; (a) its 
relative position within the range gate, and (b) its associ- 
ative conditional probability of acquisition. An ASCM 
detection/acquisition system with a multi-target discrimina- 
tion ability has a higher probability of acquiring the HVT 
than that of acquiring a LVT nearby. The model represents 
this discrimination by incorporating multiple Sweep Widths 
and conditional acquisition probabilities associated with 
the ASCM system and target set. In practice. Sweep Widths 
and conditional acquisition probabilities may be difficult 
to estimate. Possible sources are ASCM system simulations, 
observational test firings, and expert opinion. 

One final assumption is made concerning the nature 
of the offensive firing. Destruction of defensive firepower 
is the initial result of missile hits on the HVT and her 
escorts. Target sinkings occur long after the ASCM salvo 
has run its course. Hence, although hits on secondary 
targets degrade defensive firepower, total target elimina- 
tion is not likely to occur during the search and detection/ 
acquisition phase of the ASCM attack ( ie. , an HVT or LVT 
that is present for detection by the first missile in the 
salvo will be present for the last missile to detect). 

3. Modeling the Target Formation 

The target formation is modeled as a single HVT 
within a bounded, spatially homogeneous Poisson random field 
of escorts. The Poisson field completely characterizes the 
number of LVT in the field ( a Poisson random variable given 
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by the field density parameter) and the placement of LVT 
within the field (uniformly randomly distributed conditioned 
on the number). It may seem that the Poisson field is 
better suited for modeling a scattering of merchant shipping 
in a sea lane, rather than a geometrically rigid Anti-air 
Warfare (AAW) screen. In reality, screens are not rigid: 
escort ships are randomly patrolling sectors, exchanging 
stations, detaching to form Anti-Submarine Warfare (ASW) 
groups, etc. , making for a dynamic structure in contrast 
with the mental image of a rigidly symmetric AAW screen. A 
desirable feature of the Poisson field is the notion that 
the number of LVT in the field is a random variable. This 
is likely to be true if the strike planner's intelligence 
estimate of enemy force composition is uncertain and 
targeting of secondary LVT is unavailable. The average 
number of escorts accompanying the HVT and the area of the 
Poisson field can be estimated with the aid of intelligence 
reports, electronic support measures, and subjective prob- 
ability encoding, thereby furnishing a reasonably good esti- 
mate of the field density parameter. For example, a field 
of 8 escorts uniformly distributed within a circle of radius 
12 nautical miles about the HVT has a density of 0. 0177 LVT 
per square nautical mile. As will be seen, such a figure 
can be utilized to select a nearly optimum range gate 
setting. 

B. ANALYTICAL DEVELOPMENT 
1. Geometry 

A few assumptions concerning the orientation of the 
Poisson field of escorts and the effective search area of 
the missiles are required in order to obtain analytical 
results: 

• The HVT is centrally located within the escort field. 

• The escort field area is large compared to the ASCM 
search area ( ie. , the search area is wholly contained in 
the field). 
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The model's analytics reference a standard 
xy-coordinate system oriented so that the y-axis lies on the 
missile flight path. The origin represents both the esti- 
mated mean of the HVT density in residence above the 
xy-plane and the center of the Poisson field of escorts. 
The effective search area of the missiles is dimensioned by 
range gate (measured along the y-axis) and Sweep Width 
(measured along the x-axis). It is important to remember 
that the effective search area can be a group of rectangles 
with different Sweep Widths for different target classes 
when a distinction is made, in which case each rectangle of 
effective search has a common y-dimension (range gate) and a 
unique x-dimension (Sweep Width). Figure 2.1 is a sketch of 
model geometry. 

2. Single Shot Penetrator 

For illustrative simplicity, the HVT true location 
density at the time of missile search is assumed to be 
circular normal. Appendix A contains a time-dependent 
general bivariate normal HVT location density as an alterna- 
tive. 

Let the random variable Y represent the y-coordinate 
of the HVT's true location, measured from its estimated 
location ( ie. , the origin). In the language of Bayesian 
statistics, the distribution of Y is the posterior distribu- 
tion of the HVT's true y-coordinate, given a point estimate 
of its y-coordinate of y = 0. In the following discussion, 
y refers to a possible value for the random variable Y. 
(The same interpretation holds for the random variable X, 
the HVT's true x-coordinate. ) 

The single shot calculation requires the following 
input constants and probabilities: 

X = LVT escort field density. 

= Sweep Width governing LVT. 

= Sweep Width governing HVT. 
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Effective Search Area 



Figure 2. 1 MODEL GEOMETRY. 
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dj = Near edge of range gate (y coordinate). 
d 2 = Far edge of range gate (y coordinate). 

■p = Missile reliability. 

= Conditional probability of acquiring a 
LVT given LVT is detected. 

= Conditional probability of acquiring the 
HVT given HVT is detected. 

Assume that the HVT is contained in the ASCM's 
effective search area. Since detected targets are processed 
from near to far within the missile seeker sensor's range 
gate, HVT acquisition requires that no LVT has been acquired 
( ie. , selected and homed-on) in the region W (y-d ^ ) . By 

X/ 

virtue of the Poisson field assumption this is equivalent to 
experiencing a free path of areaW^Cy-dj) in a in a two- 
dimensional filtered Poisson process of rate . 

Define the following events: 

D: HVT acquisition. 

E: Filtered escort free path to the HVT. 

The probability of HVT acquisition is expressed as 
d2 W /2 

PCD) = / rP P(ElY=y) f(x,y) dxdy, (2.1) 

di -W, /2 
^ h 

where 

-Cx^+y^ ) /2o^ 

f(x,y) = e 1 , (2.2) 

2iTa^ 



is the probability density function describir- 
of the HVT, and 



P(E|Y=y) 



-XW P (y-di) 



the location 



(2.3) 



is the conditional probability of event E, given the true 
range coordinate of the HVT. In terms of these components, 
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PCD) 



rP, e 
h 



do VL/2 

/ /*’ 
d, -W. /2 



-XW^P^(y-di) -Cx^+y^)/2a2 
e 



XW P di 

rP- e ^ {$(W, /2a) - $(-W^/2a)} 

h h h 



d 2 -^W P y -y2/2a2 
X / e e dy . 



1 



/2it< 



dxdy , 
2ira^ 



The integral is evaluated by completing the square in the 
exponent and adjusting the limits of integration. 

Let K = XW P , z = yiil£_ ; 

i 1’ a ’ a ’ 



PCD) 






Kd 2 ( K cr ) ^ 

e { <l>(d 2 + Kcy^ ) 

X 

a 



^ ( d 1 -f ) I, 
a 



- 



Kd 1 C K a ) ^ 

e AO A4> ' , 

X y 



C2 .4) 



where 



b 

AO = / 4>Cu) du 

a 



is the integral of the standard normal probability density 
function evaluated at the upper and lower limits shown. 
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3. Multiple Penetrators 

The probability that at least one of m penetrating 
missiles acquires the HVT is computed by conditioning on the 
number of LVT in Wjj^(y-di) of the effective search area. The 
solution' is reached by integrating the HVT density over the 
effective search area, after first removing the condition on 
the number of LVT "in the way". 

Let N, a random variable, be the number of missiles 
that acquire the HVT. The conditional probability that at 
least one of m missiles acquires the HVT is: 



where L is the number of LVT in W^Cy-d^). 

The condition on LVT is removed using the knowledge that the 
assumed distribution of LVT escorts is Poisson with spatial 
density parameter X . Expansion of the binomial term and 
rearrangement shows that 



P(N>0lL=i, Y = y) = 1 - (1 - rP^^C l-P^, ) ^ 



(2.5) 



oo • 

P(N>0|Y = y) = 1-1 (1 - rPj^d-Pj, 

i= 0 L 



OO 



-XW^(y-di) 



(XWj^(y-di))^ 



1-1 E 

i=0 k=0 . 



m 




(-l)^(rPj^(l-P^)^)^ 



-XW^ (y-di ) 



( XW^ (y-di ) ) 



1 



i! 
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1 



m 

- Z 

k =0 




(-rP^)’' 



-XW^ (y-di ) 
e 



X I (d-P ( XW^(y-di) )^ ^ 

i =0 ^ i! 



= 1 



m /m\ ^ -XW^(y-dj)(l - (1-P^)^) 

Z Vkj (-rP, ) e 



( 2 . 6 ) 



Finally, to remove the condition on the location of the HVT, 
equation (2.6) is multiplied by the HVT location density and 
integrated using the same technique employed for the single 
shot calculation. 



Let 6 (k) = XWj^Cl - ); 



d 2 W, / 2 m 
P(N> 0 ) = 1 - / ^ 

di -W, /2 k =0 L 



il] 



(-rP, ) e 
h 



- 5 (k) (y-di ) 



- (x^+y^ ) /2o^ 

X e dxdy 

2 IT a 



m 

1 - A$ E 
^ k = 0 



f m\ . 6 ( k ) d 

ik (-rP, ) e 
h 



1 



dz - 6 (k)y -y ^ / 2 

X / e e dy , 

/2Tra 
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m 



(") 



.k 



1 - AO Z Vk/ (-rP, ) e 
k=0 



6 (k)d^+^( 6 (k)a ) 



AO ' , 

y 



m 



(") 



.k 



AO E Vk/ (-rP, ) e 
k=i 



6 (k)d]^+i 5 ( 6 (k)a ) 



AO ' , 

y 



(2.7) 



where 

AO ' = O(do+6(k)o2) - 0(di +6(k)a2) , 

y 5 

Expressions (2.4) and (2.7) allow the analyst to determine 
an optimal near edge of range gate setting, d^ , ie. one that 
maximizes HVT acquisition probability, given the various 
operational parameters of the ASCM. 

4. The Fall of Multiple Penetrators 

More information can be obtained from the condi- 
tional distribution of the fall (distribution) of multiple 
penetrators, given that the HVT is contained in the 
missile's effective search area. Fall of shot terminology 
(e.g. , dispersion) is normally associated with the firing 
assessment of ballistic projectiles, but can be adopted for 
ASCM if one can estimate the probability of "absorption" of 
an ASCM into each target class. Simply stated, a missile 
which penetrates the HVT area-defense will either have the 
HVT in its projected search area or not. If not, it is no 
longer a threat to the HVT, though there is a probability, 
usually small, that it will home on an escort. For the 
missiles that are a threat to the HVT, it is possible to 
estimate how they are distributed between HVT, escorts, and 
misses. 
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The acquisition decision of a single, conditionally 
independent penetrator is likened to a single trial of a 
multinomial experiment which, given the HVT is in the effec- 
tive search area, has as its possible outcomes: 

A: ASCM acquires a LVT on the near-side of the HVT. 

B: ASCM acquires the HVT. 

C: ASCM acquires a LVT on the far-side of the HVT. 

D: ASCM fails to acquire a target. 

The probability formulas of the above events are 
initially stated with the condition of the HVT's 
y-coordinate in the effective search area. 

Let Ad = d 2 - dj , and condition on Y=y; then 

-K (y-di ) 

P(A1 Y=y) =r(l-e ), (2.8) 

-K (y-d 1 ) 

P(B|Y=y) = e , (2.9) 

-K(y-dj) -ic(d 2 -y) 

P(C|Y = y) = r(l-P^)(e )(1 - e .) , (2.10) 

-KAd 

P(D|Y=y) = (1-r) + r(l-P^) e . (2.11) 

Derivations of the mean, variance, and covariance of the 
multinomial distribution's binomially distributed marginal 
distributions are given in many basic probability texts. 
Removing the condition of the y-coordinate of the HVT within 
the effective search area and normalizing with respect to y 
produces the following marginal expectations of the condi- 
tional multinomial events ( ie. , the expected number of 
missiles absorbed by each target class, given the HVT is in 
the effective search area). 
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Let be the random number of penetrators that 

"absorbed" into cell i, where M. + = m. The 

A D L D 

pected value of is mp^, where 



= P(A) = r 
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The general procedure for computing a marginal variance is 
shown below, followed by the equation of the variance of 
event B; the number of HVT acquisitions, given the HVT lies 
within the effective search area. In general, 

Var(Y) = E{ Var(Y|X) } + Var{ E(Y]X) }. 

Therefore , 

Var(M^) = E{mp^(l - p^(y))} + Var{mp^(y)}, 

= E{mp^(l - p^(y))} + m2|*E{p?(y)} - (E{mp^(y)}) , 

2 

= m(m - l)E{p?(y)} + mE{p^(y)} - m2(E{p^(y)}) , 

, 2 

= m(m - l)E{p?(y)} + mp^ - (mp^) . (2.16) 

Equations (2.8) through (2.9) can be substituted into equation 
(2.16) in place of p^(y) and the corresponding expressions 
(2.12) through (2.15) substituted in place of p^. The expect- 
ation of p?(y) remains to be calculated and can be accomplished 
using 

da / d2 

E{g(p.(y))} = / g(p.(y)) <j)(y) dy / / <J)(y) dy , 

di ^ / di 

the expectation of a function of p^(y) , given the HVT is in 
the missile’s effective search area (hence the normalization 
quotient ) . 

Calculation of the variance of Mg, the number of 
ASCM that acquire the HVT, given the HVT is contained in the 
missile's effective search area is as follows; 
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(mpg) , 



where 



A(j) ' ' = $( d? + 2Ka^ ) - $( di +2 ko^ ) . 



It is important that the effect of the range gate 
decision variable on HVT acquisition probability can be 
explored using the model's analytical equations. An illus- 
trative numerical example is presented in the next chapter. 



25 



III. OPTIMAL RANGE GATING 



The near edge of the missile seeker sensor's range gate, 

, is an influential decision variable. As the range gate 
increases in length, more targets are likely to be detected 
and processed for acquisition, making the ASCM's task of 
selecting the HVT more difficult. As the range gate 
narrows, the probability of capturing the HVT in the effec- 
tive search area ( ie. , detecting the HVT) decreases. If the 
far edge of the range gate, d 2 , has been chosen such that 
the subsequent choice of d^ brackets the mean of the HVT 
targeting density along the missile flight path, equation 
(2.7) can be used to solve for a value of d^ that maximizes 
the probability that at least one penetrating missile 
acquires the HVT, given m penetrators are realized. If an 
estimated distribution of the number of penetrators were 
available, the law of total probability could be used to 
obtain an unconditionally optimal or "weighted" value of d^ 
based on the conditionally optimal values of d^ obtained 
from graphing equation (2.7) for m penetrators, m = 
l,2,...,s, where s is the total salvo size. Or, in the 
situation where only an expected number of penetrators is 
available, equation (2.7) can be graphed once for the 

optimal value of d^ based on the expected value of m. The 

numerical example that follows illustrates how the optimal 
choice of d^ is affected by ASCM multi-target discrimination 
ability, the number of penetrators, and HVT targeting accu- 
racy. 

The example uses the following hypothetical input param- 
eter values: 

0 = 3.0 Nm 
X = 0. 0177 LVT/Nm2 

W, = 6 Nm 
n 
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= 6 Nm 
P-h = 0. 97 
P* =0.4 
r = 0. 99 
d 2 = +5.0 Nm 

Figure 3. 1 is a plot of a single shot penetrator's HVT 
acquisition probability versus di , where zero on the hori- 
zontal axis is the expected location of the HVT along the 
flight path. The range gate [ dj , d 2 ] = [-5.1, +5.0] maximizes 
the single shot probability of HVT acquisition. Figure 3.2 
contrasts the curve of Figure 3.1 (curve A) with a curve 
generated from identical input parameter values excepting 
conditional acquisition probabilities, which have been set 
equal to 0.97 (curve B). The relatively gradual ramp of 
curve A is representative of a "more forgiving" missile 
( ie. , one that has a better multi-target discrimination 
ability) . 

Figure 3. 3 demonstrates the increase in HVT acquisition 
probability from firing a salvo and realizing four penetra- 
tors versus firing a single shot with the hope of realizing 
one. While it is difficult to spot a noticeable difference 
in d^ associated with the maxima of each curve, the actual 
difference is about 1. 7 nautical mile, with the four pene- 
trators having the "longer" optimal range gate. The result 
is a consequence of maximizing the probability that at least 
one missile acquires the HVT. The marginal "single shot" 
probability of HVT acquisition of each missile is effec- 
tively reduced while the chance at least one acquires the 
HVT is increased. This phenomenon becomes more apparent in 
Figure 3. 4, ten penetrators versus a single shot penetrator. 

Figure 3.5 is a plot of number of penetrators versus HVT 
acquisition probability where the conditionally optimal 
range gate has been used at each increment of the indepen- 
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dent variable. The marginal return of the first few pene- 
trators is striking. Figures 3.7 and 3.8 demonstrate the 
marked effect that an increasing HVT location error sigma 
has on HVT acquisition probability in both the single shot 
penetrator and salvo/multiple penetrator cases. As sigma 
increases, it becomes apparent that the simple firing 
strategy may very well be out-performed by a strategy that 
involves firing a spread of missile salvos ( assets permit- 
ting) , an interesting hypothesis for future research. 

Table I summarizes conditional fall of shot expectation 
and standard deviation for a single shot penetrator and 
various numbers of multiple penetrators over selected values 
of sigma. Again, the fall of shot for events A, B, C, and D 
is conditioned on the HVT's containment and location within 
the missile's effective search area. As sigma increases, 
the optimal range gate opens, the probability of containing 
the HVT in the effective search area decreases, and, in the 
case of salvo/multiple penetrators, the conditional fall of 
shot expectation becomes more evenly distributed between the 
HVT and the LVT on the near -side. It should be noted that 
the total expected threat to the escorts cannot be deter- 
mined by simply adding the expectation of events A and C, 
because there is a chance that some escorts may lie in the 
missile's effective search area though the HVT does not. 
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SINGLE SHOT PENETRATOR 




Figure 3.1 SINGLE SHOT P( HVT ACQ. ) VS. RANGE GATE. 
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range: gate: 



SINGLE SHOT PENETRATOR 




(NOlilSinODV lAH)d 



Figure 3.2 THE EFFECT OF MULTI-TARGET DISCRIMINATION. 
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RANGE GATE 



SINGLC SHOT PENETTRATOR VS. FOUR PENETRATORS 




Figure 3.3 ONE VS. MANY. 
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RANGE gate: 



SINGLE SHOT PENETRATOR VS. TEN PENETRATORS 




Figure 3.4 RELAXING THE RANGE GATE DECISION VARIABLE. 
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RANGE GATE 



NO. PENETRATORS VS. HVT ACQUISITION PROBABILITY 




Figure 3.5 INCREASING PENETRATORS. 
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PENETRATORS 



SINGLE SHOT PENETRATOR 




(NOlilSinOOy lAH)d 



Figure 3.6 TARGETING ACCURACY AND THE SINGLE SHOT. 
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RANGE GATE 



TEN PENETRATORS 




Figure 3.7 TARGETING ACCURACY AND MULTIPLE PENETRATORS. 
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RANGE CATE 



TABLE I 

CONDITIONAL FALL OF SHOT 



SINGLE SHOT PENETRATOR 



a 


(di ,d2 ) 


P( HVT)* 


Expectation ( Std. dev. ) 
ABC 


of events; 
D 


1. 5 


-3. 0, +5. 0 


0. 932 


0. 119 


0. 844 0. 004 

( 0. 362) 


0. 031 


3. 0 


-5. l,+5. 0 


0. 619 


0. 188 


0.777 0.004 

( 0. 415) 


0. 029 


9. 0 


-12. l,+5. 0 


0. 162 


0. 315 


0.654 0.005 

( 0. 475) 


0. 024 


FOUR PENETRATORS 










a 


( , d2 ) 


P( HVT)* 


Expectation ( Std. dev. ) 
ABC 


of events 
D 


1. 5 


-4. 0,+5. 0 


0. 950 


0. 614 


3.24 0.019 

( 0. 801) 


0. 121 


3. 0 


-6. 8, +5. 0 


0. 642 


0. 947 


2.92 0.018 

( 1. 000) 


0. Ill 


9. 0 


-16. 0,+5. 0 


0. 175 


1. 57 


2.31 0.022 

( 1. 310) 


0. 088 


TEN PENETRATORS 










a 


(d^ ,d^ ) 


P( HVT)* 


Expectation ( Std. dev. ) 
ABC 


of events 
D 


1. 5 


-4. 8, +5. 0 


0. 953 


1. 81 


7.84 0.046 

( 1. 38) 


0. 295 


3. 0 


-8. 3 , +5. 0 


0. 648 


2. 81 


6.87 0.043 

(1.89) 


0. 268 


9. 0 


-19. l,+5. 0 


0. 181 


4. 54 


5.19 0.053 

(2. 66) 


0. 206 



* P(HVT) = "The prob. of HVT containment in the search area". 
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IV. SUMMARY 



The number of ASCM area defense penetrators necessary to 
achieve a required minimum number of HVT acquisitions is a 
random variable. Tables or nomographs containing number of 
penetrators, optimal range gate, expected number of HVT 
acquisitions, and standard deviation of HVT acquisitions 
versus escort field density could be generated for a set of 
different posterior distributions describing HVT location 
over a range of Poisson field densities. It would be neces- 
sary to accompany these "conditional" estimates with a prob- 
ability statement of HVT containment in the effective search 
area. The model in its present form serves best as an 
example of how a probability model incorporating a rela- 
tively simple analytic form can be used for OTH strike plan- 
ning. Complex firing situations and target coverage 
problems involving multiple launches from a consortium of 
strike groups could be addressed by expanding the simple 
analytical framework of the model. The equations of a well- 
tailored (and well-validated) follow-on model could be 
contained in an interactive computer program which accepts 
periodic targeting updates and in turn furnishes the strike 
planner with new estimates of salvo size and optimal range 
gating to achieve a desired probability of HVT acquisition 
based on the number of area defense penetrators. 
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APPENDIX A 

UPDATING THE HVT LOCATION DENSITY 



The time delay measured from the generation of the HVT 
location density to the time the missile commences its 
search will often allow considerable movement of the enemy 
formation. A postulated bivariate normal target velocity- 
time vector can be used to revise the "old" HVT density. 
The approach does not win "something for nothing" since an 
assumption concerning enemy intentions is required. 
Constant target motion is a reasonable assumption if (a) the 
enemy has not been alerted by the opposition's tracking/ 
targeting methods and (b) timely early-warning of the 
impending strike is unlikely. A Brownian motion or fleeing 
datum assumption would be more appropriate in the alerted 
enemy situation. Dispersion of the enemy formation is also 
a possibility in the event of a probable nuclear exchange. 
The constant motion model developed below may prove valid 
when offensive targeting is reasonably passive and covert. 

Let OX, the HVT's vector of movement to its future posi- 
tion measured from the mean of the "old" HVT location 
density, be the resultant of the sum of two independent 
bivariate normal vectors: a position error, OX' , and a 

velocity-time vector X'X, (see Figure A. 1 ) . 

The distribution of OX' is given by the HVT density, 

OX' Uy 5 XY ^ 

while X'X is distributed as follows, 

X'X --BVNCiiyt , Vyt, agt2, ajt2, 
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X 




ox = OX' + X'X 



Figure A. 1 THE SUM OF TWO INDEPENDENT RANDOM VECTORS. 
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The distribution of OX is gleaned from its moment gener- 
ating function formed by the product of the individual 
moment generating functions of OX' and X'X. The moment 
generating function of the bivariate normal distribution is: 

= exp{siy^ + S 2 yy + 2 pSiS 20 yaY 

+ sloy^ } • 

exp{siVy + S2yy + hisja^ + 2pyySiS20yay 
+ s^o^)}, 

exp{siPyt + S 2 yyt + 

+ 2p^ySiS2aytOyt + s^ojt^)}, 

OX ' ^ ^ 1 ^ 2 ^ 5 

= exp{Si(Uy + Vut) + S2(Py + Viyf) 

+ m S^COy + O^t^) + 2 s S 2 ( PyyOyOy 
^ + ajt^) )}. 

By inspection, 

OX'^BVNCyy + y^t, yy + y^t, aj + ogt^ , aj + ajt^ , 

‘^XY^X'^Y 



Therefore, 

^OX ' ^ ® 1 ’ ^2 ^ “ 



^X *X^^^ ’® 2 ^ 



and, 

m^^Csi ,S 2 ) = m 
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The distribution of OX can be used to adjust the launch 
bearing of the ASCM salvo and update HVT postional uncer- 
tainty to the estimated time of missile search. Model geom- 
etry becomes more complicated when a general bivariate 
normal HVT location density is incorporated instead of a 
circular normal one since a unique orientation of the effec- 
tive missile search area results. 
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